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Abstract 

The dimensionally regularized massless non-planar double box Feyn- 
man diagram with powers of propagators equal to one, one leg off the 
mass shell, i.e. with p\ = ^ 0, and three legs on shell, = 0, i = 
2,3,4, is analytically calculated for general values of and the Man- 
delstam variables s, t and u (not necessarily restricted by the physical 
condition s + t + u = q^). An explicit result is expressed through 
(generalized) polylogarithms, up to the fourth order, dependent on ra- 
tional combinations of q'^,s,t and u, and simple finite two- and three 
fold Mellin-Barnes integrals of products of gamma functions which are 
easily numerically evaluated for arbitrary non-zero values of the argu- 
ments. 
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1. Feynman diagrams with four external lines contribute to many important phys- 
ical quantities. They are rather comphcated mathematical objects because depend 
at least on two independent Mandelstam variables. In the pure massless case with all 
end-points on shell, i.e., pf = 0, i = 1, 2, 3, 4, the problem of the analytical evaluation 
of two-loop four point diagrams (when the planar and non-planar double boxes shown 
in Fig. 1 are mostly complicated), in expansion in e = (4 — (i)/2 in the framework 




Figure 1: The planar and non-planar double boxes 

of dimensional regularization with the space-time dimension as a regularization 
parameter, was completely solved during last year in 0, |, ||, ^. The corresponding 
analytical algorithms were quite recently applied to the evaluation of two-loop virtual 
corrections to e~^e~ — > fi~^fi~ and Bhabba scattering and quark scattering [0]. 

In the massless case with one-leg off shell, pf = q'^ 7^ 0, p'^ = 0, i = 2,3,4, 
which is relevant to the process e~^e~ — >• 3jets (see, e.g., 0), the planar double box 
diagram with powers of propagators equal to one has been analytically calculated in 
for general values of and the Mandelstam variables s and t. An explicit result 
is expressed through (generalized) polylogarithms, up to the fourth order, dependent 
on rational combinations of g^, s and t, and a one-dimensional integral with a simple 
integrand consisting of logarithms and dilogarithms. This result represented through 
so-called two-dimensional harmonical polylogarithms was confirmed in |]T0| where the 



method of differential equations |Tl], was applied and some other planar master 
Feynman diagrams with one leg off shell were also evaluated. 

There are two different non-planar Feynman diagrams (see Fig. 1) with all powers 
of propagators equal to one: when 7^ (or p^ 7^ 0) and when p| 7^ (or P4 7^ 0). 
The purpose of this paper is to analytically evaluate the first of these two master 
integrals. As in the pure on-shell case (see [Q) it turns out that it is natural to 
consider non-planar double boxes as functions of s,t,u = {pi + PiY and q^ = pi 
not necessarily restricted by the physical condition s + t + u = q"^ which does not 
simplify the result. We shall present an explicit result for the considered master non- 
planar double box in terms of (generalized) polylogarithms, up to the fourth order, 
dependent on rational combinations of q'^,s,t and u, and simple finite two- and three 
fold Mellin-Barnes (MB) integrals of products of gamma functions which are easily 
numerically evaluated for arbitrary non-zero values of the arguments. 

To arrive at this result we straightforwardly apply the method of refs. we 
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start from the alpha-representation of the double box and, after expanding some of 
the involved functions in MB integrals, arrive at a five-fold MB integral representa- 
tion with gamma functions in the integrand. Then we use a standard procedure of 
taking residues and shifting contours^ to resolve the structure of singularities in the 
parameter of dimensional regularization, e. This leads to the appearance of multiple 
terms where Laurent expansion in e becomes possible and provides our result. 

Finally we discuss the problem of systematical evaluation of general double boxes 
with one leg off shell and the most general class of functions that can appear in the 
results. 

2. The alpha representation of the non-planar double box with one leg off shell, 
Pi = 9'^ 7^ differs from its on-shell variant (see, e.g., eqs. (4)-(6) of ref. 0]) by a 
term depending on inserted into one of the functions involved: 

F(s,t,u,g2;e) = -r(3 + 2e) (m'^/^^ / den.../ da^5 {^ai - l) D^+^' A'"^-^' , 

(1) 

where 

D = («! + ^2 + «7)(«3 + «4 + «5 + «6) + («3 + «5)(«4 + «6) , (2) 

A = [aia2{a3 + a4 + + a^) + aia^a^ + a2a3aQ]{—s) 

+a5aQaj{—t) + asa^aj^—u) + ajla^a^ + 02(03 + 0:4 + "5 + Ci6)]{—q^) . (3) 

By introducing six times MB representation 

= — / dw —T(iy + w)T(-w), (4) 

in a suitable way, one can take all the parametrical integrals in gamma functions. 
Fortunately, one of the MB integrations is then explicitly performed by use of the 
first Barnes lemma and we arrive at the following 5-fold MB integral: 

(m^/2)'r(-e)2 1 r (q2y 

^^^'^'^'^^ = - r(-l - 3e)r(-2e)(-.)3+^^ i2.rf I d^d^id^^d.,dz2 (^-j 




ty^ /uy r(3 + 2e + t; + Wl + W2 + Zi)r(3 + 2e + wi + W2 + zi + Z2) 
s) \s) T{'i + 2e + wi+W2 + zi)T{2 + wi + W2 + zi + Z2Y 



xr(2 + e + wi + W2 + zi + Z2)T{1 + v + wi + W2)T{1 + wi + zi)V{l + W2 + zi) 

Xr(l +Wi + Z2)T{1 +W2 + Z2)T{-2 -2t-v ~Wi~W2- Zi) 

xr(-2 -2e-wi-W2- Z2)T{-v)V{-w{)V{-W2)T{-zi)V{-Z2) . (5) 

It differs from its analog for = (see eqs. (8)-(10) of ref. [Q) by the additional 
integration in v. This variable enters only four gamma functions in the integrand. 

^The procedures presented in and |4| are a little bit different. In the latter variant, one 
systematically uses integration in MB integrals in straight lines along imaginary axes. We prefer to 
use the former variant. 
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When taking minus residue of the integrand in v one reproduces the corresponding 
4-fold MB integral in 0]. Integral @ is evaluated in expansion in e, up to a finite 
part, by resolving singularities in e in the same style as in 0. We obtain and confirm 
as a by-product the on-shell result of ref. [Q. 

To present the final result let us turn to the variables x = s/q"^, y = tjq^ and 
z = u/q^\ 

We obtain 

1 3 

f4{x,y,z) = -{l-x-yz)--{z + y) , (7) 

Ux,y,z) = -^[3{l + x + y + z) + {2-2x-3y-3z-2yz% 

+ {2 - 2x + y - 3z + 2yz)ly + {2-2x-3y + z + 2yz)l,] , (8) 
f2{x,y,z) = 3ix-l)U2ix) + 2{l-y)zU2iy) + 2yil-z)U2iz) 



1 r 



{2 - 2x - 3y - 3z - 2yz)ll + {2 - 2x + y - 3z + 2yz)l 



+ 4 

+ {2~2x~3y + z + 2yz)ll\ - 3(1 - x)U^ + 2z{\ - y)lyly + 2y{l - z)lj, 

+ ^ [{2 - 2x + y - 3z + 2yz)ljy + {2 - 2x + z - 3y + 2yz)lJ^ 

3 

+ {2-2x + y + z- 2yz)lyh] + -[{I - x + y + z% + {I + x - y + z)ly 

+ {l + x + y- z)U] + —(30(1 -x) + ll{y + z) + 2Qyz) + 3{l + x + y + z) , (9) 

48 

h{x,y,z) = 3{l - x){2U^ {x) + bU^ {I -x)) + 2{2-2x + Ay -z-3yz)U^{y) 

-3(1 - y)zU^ {l-y) + 2{2-2x-y + Az- 3yz)U,, {z) - 3y{l - z)U^ (1 - z) 

-3(1 - x)(5L - 2(1 + ly + h))U2 {x) - 15(1 - x%U2 (1 - x) 

+ ((1 - y)(6 - Azl^ + 3zTy + Azl,) -(2-2x + Ay + z- ^yz)ly)U2 [y) 

+ ((1 - 2) (6 - Ayh + 3yl, + 4^/^,) -{2-2x + y + Az- byz)l,)U2 {z) 

+3(1 - y)zlyU2 {l-y) + 3y{l - 2)LLi2 (1 - z) 

-^(2 - 2x - 3y - 3^ - 2yz)ll + 3(1 - x)^/^ - 15(1 - x)?J^ 
-}-{2-2x + y-3z + 2yz)ll - 2(1 - y)zlyll + 3(1 - y)zllly 
-^{2-2x-3y + z + 2yz)ll - 2y{l - z%ll + 3y{l - z)Fj, 



1 

~2 

+ {2-2x + y + z- 2yz)lyh{ly + -{2-2x + y + z- 2yz)ljyl 



[{2-2x + y-3z + 2yz)ljy{l^ + ly) + {2 ~ 2x - 3y + z + 2yz)lJ,{h + h) 
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+ 6(1 - X)lJ^{ly + h) + 4(1 - y)zlyly{h " I x) " 4?/(l - Z)hh{l^ - ly) 

3 r 

--[{l-x + y + z)ll - 4(1 - x%U + {l + x-y + z)ll - 4(1 - y)lyly 
+ {l + x + y- z)ll - 4(1 - z)Uh 

-3((1 - x-y + z)ljy + {1 + x - y - z)lyU + {I - x + y - z)lJz) 

-— [(30 - 30x + lly + 11^ + 2Qyz% + (30 - 30x + 3?/ + ll;z - lAyz)ly 

-60(1 - x% +12(1 - y)zly + (30 - 30x + lly + 3^ - Uyz)U + 12y(l - z)l, 
-Q[{l-x + y + z% + {l + x-y + z)ly + {I + x + y - z)h] - 12(1 + x + y + z) 



03) 
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[52(1 - x) - 63(y + -2) - 40?/z] (11 - 5(a; + y + z)) 



(10) 



where = In a and = ln(l — a) for a = x,y, z. Moreover Li^ (z) is the polylogarithm 
|]13| and (in the next formula) 



ly+b-i ri \rf -\t) \n\l - zt) 



{a-l)\b\ Jo 



dt 



(11) 



the generalized polylogarithm |]14| . 

The finite part /o = /o + /o consists of two pieces: 

/o(x, ^) = 30(1 - x)S2,2{x) - 2(1 -y){6- 7z)S2,2{y) " 2(1 - ^)(6 - 7y)S2,2{z) 
+6(x — yz)Li4 iyz/x) — 12(1 — x)Li4 (x) 

-2(3 - 6x - 3y - z + 4y2)Li4 (y) - 2(3 - 3x - y - 3z + 4:yz)Ui {z) 
+69(1 - x)Li4 (1 - x) - 3(1 - ?/)zLi4 (1 - y) - 3?/(l - z)Li4 (1 - z) 

-i [(14 + lOy -7z + yz)U2 {yf + (14 -2x-7y + lOz + i/z)Li2 (2)'" 

-6(1 - x) [(2 - 54 + 2/y + 2/,)Li3 (x) + (4 + 5/^, + 5/,)Li3 (1 - x)" 

- [4(2 + x + Ay~z- 3yz% + 2(1 + 5y - 5yz)ly + 4{2 - x + 4y + z + 3yz)l^ 
+2(1 -y){Q- 7z)ly + 12(1 - y)] Lig (y) 

-(6(1 - i/)(4 - + hz) + (20 + 4y - 17z + lli/;2)/j,)Li3 (1 - y) 

- [4(2 + X - y + 4^ - 3t/z)/^ + 4(2 - 2x + ?/ + 4^ + 3yz)/j/ + 2(1 - x + 5z - 5?/z)/^ 
+2(6 - 7y)(l - ^)L + 12(1 - z)] Lis (^) - [6(4 - l^y + /,y)(l - z) 

+ (20 - 2x - 17y + 4z + llyz)Q Lis (1 - -2) + Q{lx - ly - lz){x - yz)U3 (yz/x) 



+ 



;(1 - X) [4{ly + h + I f + TT^ + 12] Lis (X) 

4(1 - y)z{k - Izf + {2-2x + Ay + z- hyz)ll + 2{2 + x + Ay + z- 5yz)U 



+2(2 -x + Ay- z + 5yz)lyl^ - lyly{lA + lOy -7z + yz) - 12(1 - y){l^ + 
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2 2 " 

+y (22|/ + Ibz - 2lyz) - 24(1 - 2/) + y Lis (|/) 

+ [4y(l - z){l:c - lyf + (2 - 2a: + y + 4z - 5yz)Z^ + 2(2 + x + y + 4z - 5yz)Za;^z 
+2(2 - 2x - y + 4z + 5yz)y^ - (14 -2x-ly + IQz + yz)l^l^ - 12(1 - z)(Z^ + ly) 

Li2 (^) 



- — (2x - 15|/ - 22;^ + 2lyz) - 24(1 - ^) + y 



+(x - yz)(3(Z^ - - l,f + 47r')Li2 (y^/x) 

+^(2 - 2x - 3y - 3z + yz)Z^ - 2(1 - x)l% + y (1 - a:)/^/! 

+1(2 - 2x + y - 3^ + 5y^)Zj + ^(1 - y)zllly - (10 + 2y - 7^ + Ayz)ll?y 

+1(2 - 2x - 3y + z + 5y;2)/^ + ^y(l - z)l% - (10 - x - 7y + 2^ + 4y^)/2/^ 

+^(2 - 2x + y - 3z - y^)/./j,(/^ + i^) + 1(2 - 2x + y - 3z + Sy^)^^^^ 

+1(2 _ 2a; - 3y + 2; - yz)lMll + Z^) + 1(2 - 2a; - 3y + 2; + Sy^:)/^/^ 

+1(2 - 2a; + y + ^ + yz)lyh{ll + Z^) + 1(2 - 2a; + y + z + 9yz)Z^Z^ 

-6(1 - x)lJx{ly + ^2)(^a; + + ^2) 

+4(1 - y)zlyly{l^ - + ly- h) + 4y(l - z)l^{l^ - ly)U{lx -ly + h) 
+ [(2 - 2a; + y + z + yz% + (2-2a; + y + ^- f)yz){ly + Q] IJyl, 
-{x - yz)% -ly- Izf^O- - y^/^) 

+(1 - a; + y + z)ll - 6(1 - x)l% + (1 + a; - y + z)ll - 6(1 - y)llTy 

+(1 + a; + y - z)ll - 6(1 - z)llU + 3((1 - x - y + z%ly{l^ + Z^/) 

+(1 -x + y - z)lJz{l^ + 1^) + {1 + X - y - z)lylz{ly + Z^)) + 6(1 - a; - y - z%lylz 

-12 [lX{ly + Q{1 -X)+ lylyik + ) ( 1 " y) + (l " z) + ) 

+7r' {^(30 - 30a; + Uz + y(ll + 74z))ll - y (1 - x%k 

+^(30 - 30a; + 3y + Uz + 98yz)/2 + 1(14 + y(10 - 7z) + z)lyly 

+ ^(30 - 30x + lly + 3^ + 98yz)ll + 1(14 - 2a; + y + 10^ - 7y^)/X 

+ 1 [(30(1 - a;) + 3y + 11^ - 62yz)ljy + (30(1 - x) + lly + 3z - 62yz)U, 
+ (30(1 - x) + 3y + 3^ + 50y^)yj - 4(/, - Zj, - /,)(a; - yz) ln(l - y^/a;)} 



+6 



(1 - X + y + z)ll - 4(1 - x%k + {l + x-y + z)ll - 4(1 - y)lyly 



+ {l + x + y-z)ll-A{l-z%l, 
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+12((1 - X - y + z)ljy + x + y - z%U + {I + x - y - z)lyl^) 

+C(3) [-^(38(1 -x) + 63(y + z)+ A{)yz% - 30(1 - x% 
4 

+-(28 - 13x + 1% - 2Az + 10?/z)/y + 2(1 - y)(6 - r^)/^, 
o 

+|(56 - 29x - 48?/ + 28z + 2Qyz)U + 2(6 - ly){l - 2)!, 

+Y [(11(1 -x)- h{y + ^))/, + (11(1 -y)- 5(x + z))ly + (11(1 - - 5(x + y))h] 
+24((1 - X + y + z)/^ + {l + x-y + z)ly + {l + x + y- z)l^) 

+48(1 + X + y + z) + —{II - h{x + y + z)) + 4(17 - x - y - -2)C(3) 

+ ^(894(1 - x) + 395(y + + 10034y^) (12) 

and the following two- and three-fold MB integrals: 

/o(x, y, z) = -7-^ / dvdw^dw2 x-i— 2^i+-izi+-2r(3 + v + + W2) 

xr{i + v + wi + W2)r(i + wi)r{-wi)r{-i - wi)r(-i - - wi) 
xr(i + W2)r(-^2)r(-i - ^2)r(-i -V-W2) 

+ 77^ I ^vdw x~'^{y^^^z-^ + z^+^y-"") 
xr(l + v)T{v)T{~v)T{l + wfT{w)T{-wfT{-l - v - w) 
+ — ^ / dt;du; x-''-'"{y^+'" + zi+'")^^^^r(l + u;)r(-l - wyT{l + v + w) 
xT{v + w)T{-v -w)[l+v + w + v{l + w)(^(l -v)~ ^(1 + w))] 

+— ^ / dwdu; [(1 + t; + w)w{{l -v){v + w)x + - 5(1 - v)vz] 

[2m) J 

X ,^^^^} r(w)r(-w)r(-i - w)t{i + v + w)t{v + w)r(-i -v-w) . (i3) 

(1 — f )f 

The integration contours in the above 3-fold MB integral are in straight lines along 
imaginary axes with —1 <Ref , Rewi,2, Re {v + "^1,2), Re {v + wi + W2) < and, in 
2-fold MB integrals, with —1 <Ref , Rew, Re {y + w) < 0. 

3. Our result (^-(0) is in agreement with the leading power behaviour when 
— > which is obtained by use of the strategy of regions of [^. It also agrees 
with results based on numerical integration in the space of alpha parameters [|l^ 
(where the 1% accuracy for the 1/e and e° parts is guaranteed). The result has 
been also checked with the value of the double box at s = which was evaluated 
by the same technique based on MB integrals. Note that F(g^,g^,t,u) is expressed 
not only through (generalized) polylogarithms, as the on-shell dobule box in 0], 
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but also through a two-fold MB integral that leads to a one-parametrical integral of 
dilogarithms and logarithms. 



The two-fold MB integrals in (13) can be converted into a form similar to the 
two- fold MB integrals present in the planar double box or, expressed through 



so-called two-dimensional harmonic polylogarithms [|T0| which generalize harmonic 



polylogarithms introduced in [|T^ . The three-fold MB integral in ([13D can be written, 



up to some two-fold MB integrals, through the following three-fold series 

where A„ „^ „2 is a linear combination of {ipin + ni^2 + 1)^ + 'ip'in + ni^2 + 2)), 
ipi^n + ni^2 + 2), and ip{ni^2 + !)• If we replace An,nx,ni by one the corresponding 
series can be summed up with a result in terms of polylogarithms, up to Lis, depend- 
ing on certain rational combinations of x,y,z. The presence of An,ni,n2 makes the 
situation much more complicated. The series can be summed up in the form of a 
two-parametrical integral of a cumbersome expression with polylogarithms. It is yet 
unclear where one of these parametrical integrations can be explicitly performed. We 
therefore prefer to leave the 3-fold MB integral in the result. This form can be easily 
used for numerical evaluation for arbitrary non-zero values of x,y and z. 

Keeping in mind the characteristic feature of the non-planar diagrams mentioned 
above one can conclude that the presence of four scales (or, three dimensionless 
variables) in the problem shows that the two-dimensional harmonic polylogarithms 
(introduced in [|l^ just for the problem of the evaluation of the double boxes) could be 
insufficient to express all the results for the double boxes with one leg off shell, and one 
can imagine the necessity of introducing three-dimensional harmonic polylogarithms. 
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